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[ Abstract 



C/3 ^ It is shown that in the case of free electron in a spatially periodic magnetic field the concept of magnetic 

translations operators is still valid and, moreover, these operators can be defined in the same way as for a Bloch 
electron in a uniform magnetic field. The results can be a useful tool in investigation of lately observed phenomena 
in 2DEG with spatially modulated density. 
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^ ■ 1 Introduction 

o; 

^ O^ , Recently there is a big interest in properties of charged particles in different inhomogeneous magnetic field profiles, 
what has been studied both experimentally ||l| and theoretically ^, ||] . Properties of the strongly correlated phases 
T— I , can be obtained from the photoluminescence spectra so much efforts are devoted to this technique 0. This 
^ problem is strongly related to Anderson localization of electrons and, therefore, is important in investigations of 
high-Tc superconductors and composite-fermions in the quantum Hall effect. Some recent experiments have shown 
^y-j ■ that variation of density in two-dimensional electron systems lead to a fictitious periodic magnetic field [|| and a 
, geometric resonance of the classical cyclotron orbit and the field period. A theory of this effect has been lately 
T-H ' proposed Zimbovskaya and Birman Q . 

0^ , The aim of this paper is to consider a possibility of introducing magnetic translation operators in the case of 

a spatially periodic magnetic field and, if this introduction is possible, to determine their form. Some general 
considerations on two-dimensional quantum systems in a singular vector potential A were presented by Aral j7| , who 
generalized some concepts to nonuniform magnetic fields for continuous quantum systems and presented a method 
of reduction to lattice quantum systems. A special class of introduced operators leave invariant the Hilbert space 
L^{TZ^) (or 1^{Z^) for lattice systems) invariant, so it is possible to construct a continuous version of a Hamiltonian 
^ ', of the Hofstader type |^ . The problem analyzed in the present work is simpler, however it is directly related with the 
O I above mentioned phenomena. It is shown that magnetic translation operators introduced independently by Fishbeck 
^ , Brown ]lO| , and Zak ||ll| for an electron in a uniform magnetic field and a periodic potential can be also applied 
in the case of a spatially periodic magnetic field. In the simplest case a free electron is considered. 

In the next section a brief summary of the main results for a uniform magnetic field is given, whereas the problem 
^ • is solved in Sec. 3. Some remarks are presented in the last section. 

2 Bloch electron in homogeneous magnetic field 

For the sake of simplicity it is assumed that electron can move in the two-dimensional plane xy, whereas the magnetic 
field H = Hz is perpendicular to it. Therefore, the components and Ay of a vector potential A are only relevant. 
In both cases, i.e. for homogeneous and inhomogeneous magnetic field, the Hamiltonian is given by the well-known 
formula 

n = ±ip-lAr + Vir). (1) 
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The only requirement is that ^(r) is a periodic function of r G 7^^. Positions of crystal nodes are determined by 
vectors R of the two-dimensional translation group T isomorphic to Z x so they will be often replaced by pairs 
{Ux, riy) of integers, i.e. by their coordinates in the crystal base {ai, ^2}. Since a point group symmetry is not taken 
into account, the square lattice with ai _L a2 and ai = 02 = a is considered. 

A form of the vector potential For a constant and uniform magnetic field the vector potential A can be written 
as linear function of the position vector r |p^ . Moreover, the coordinates and Ay can be chosen in such a way 
that they do not depend on x and y, respectively.^ So in the most general case one obtains 

A'x = oty and Ay = [3x . (2) 

Real numbers a, (3 have to satisfy 

for a given magnitude H of the magnetic field. Of course, this form includes the antisymmetric gauge (r x H)/2 for 
a = — /3 — ~H /2 and the Landau gauge for (3 — H^a — 0. 

The symmetry of a problem The periodicity of V{y) yields the two-dimensional translation group T ~ Z x Z 
to be the symmetry group of the Hamiltonian. Therefore, one can introduce concepts of quantized quasi-momentum, 
energy bands and Bloch electrons. Irreducible representations of T label energy levels of Ti. However, in the 
presence of the magnetic field Hi, at least except for special values of H, projective representations have to be used 
and their factor systems depend on the magnetic flux through the unit crystal cell. The periodic boundary conditions 
imposed on (projective) representation T(R) leads to magnetic flux quantization and the concept of magnetic cells (or 
magnetic periodicity) [|lO[ ^ . The periodic boundary conditions are also responsible for choosing the magnetic field 
perpendicular to a crystal plane [ pO| , pT| , |T^. It should be underlined that in Zak's approach projective representations 
of T are replaced by vector representations of a central extension of T by a group of factors included in U(l). Of 
course, these two approaches are equivalent. 

Projective representations Operators of a projective irreducible representation of T which commute with TL 
can be chosen as §, |ll|, |l| 

r(R) = exp[R • (p - ^A')] , (3) 

where A' is a vector potential associated with A and, in this simplified considerations, given as 

A'^ = I3y and A'y = ax . (4) 

The commutation T(R) with Ti. follows from the fact that coordinates of tt = p — |A commute with coordinates of 
tt' = p — I A' Moreover, the commutators [tt^;, tTj,], [tt^, tt^] etc. are numbers, so they commute with any 

other operator. This fact is very important in the derivation of a factor system for the representation r(R). This 
factor system and the group-theoretical commutator r(R)r(R')T(R)~^r(R')~^ depend only on the magnitude iJ, 
not on a form of the vector potential. However, working with a local gauge A^ one can change a factor system, but 
the commutator is unaffected (is gauge-independent) [ p^ . 

Movement of a Bloch electron In the case of homogeneous magnetic field an electron moves around the cyclotron 
orbit with coordinates of center given by operators 7r(^ , tt^ . This movement is quantized and is related to the broadening 
of the Landau levels p^ . 

3 Inhomogeneous magnetic filed 

In order to keep symmetry described the translation group one has to assume that non-constant terms in the magnetic 
field magnitude are periodic with respect to x and y coordinates. Therefore, tt^, the generalized kinetic term of the 
Hamiltonian (|l|), is invariant under R e T and, at least in the first-order approximation, the potential V'(r) can 

^This condition is stronger than the radiation gauge V ■ A = 0. 
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be omitted. Applying the Fourier transform a periodic magnetic field can be written as a sum the sine and cosine 
functions sin k^x, sin kyy (or exp(ik • r), if one prefers the complex analysis), so, to simplify a problem, it is assumed 
that 

H = Ho + Hi [cos(A;a;) + cos{ky)] , (5) 

where k = kx — ky — 27r/a. 



3.1 The vector potential 

For the magnetic field given by ^ the vector potential A can be chosen in many gauge invariant forms. For the 
sake of clarity the simplest form in the radiation gauge (V • A = 0) is assumed 

Ax = ay ~ — sin{ky) , 
Hi 

Ay — f3x -\- —— s\n{kx) , (6) 

K 

with (3 — a = Ha. 



3.2 The symmetry group 

Due to the periodicity of H the symmetry group is still T ^ Z x Z. Its irreducible representations should label 
eigensapces of TL. There is, however, one problem: 'Does it suffice to consider projective representations or, due to 
higher terms of x and y in A^^Ay, and H, some more complex structure should be used?' Within the frame of Zak's 
approach it may mean that one has to investigate non-Abelian extensions, for example. It seems that when for A 
being linear function of r the second cohomology group (related with projective representations or, equivalently, with 
central extensions) comes into play, then for A^, Ay being second order functions of x,y the third cohomology group 
should be considered etc. There are some hints from mathematics and physics which indicate that it suffices to limit 
to projective representations (i.e. the second cohomology group). 

At first, we are interested in A and H expressed by the (co-)sine or other periodic functions with an infinite 
Taylor expansion, so — if the order of a cohomology group depends of the order of functions — cohomology group of 
the infinite order should be taken into account. Moreover, factor systems of projective representations in the case 
of a homogeneous magnetic field depend on the magnetic flux through the unit crystal cell. This quantity is always 
a number calculated as an integral over the unit cell of the product H{x, y)dxdy (or as an integral over the edges 
of a cell /A(r) ■ dr §). Such an integral depends only on the constant term Hq since the integration of periodic 
termsQ gives zero. The cohomology group of order n demands considerations of the nth co-boundaries and co-cycles, 
which involve n lattice vectors. It seems it would be necessary in the case of a hypothetical 'field' being a tensor 
of rank n {F^yz^Gxyzt, etc.) — in the case of magnetic field, which can be described by the tensor Hxy, the most 
important are loops (drawn in two dimensions) encircling the magnetic flux (a one-dimensional object). At last, 
there are no well-investigated group-theoretical generalizations of non-Abelian extensions corresponding to fourth, 
fifth etc. cohomology groups (at least the author is not aware of such considerations). 



3.3 Projective representations 

In this section a form projective representation T of T is derived. This representation should commute with Ti and 
the results for the uniform magnetic field have to be revealed in the limit Hi —>■ 0. 
Let us assume that operators of a projective representation have the following form 

T(R) = exp(^-i7r'.R^ , (7) 

where a form the operator tt' depends on A and in the limit _ff i we have tz' = p — (e/c)A' with A' defined in 
the previous section by Eq. (^). Due to the periodicity of V{r) (moreover, it is assumed V{r) = in the simplest 
approximation) it is enough to calculate a commutator [T'(R),7r^]. If tTx and tTj, commute with T'(R) then also H. 

^In this place a 'periodic ierm ' denotes all but a constant term of the Fourier transform. 
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does. The condition [T(R),7ra:] = [T(R),7rj,] = allows also labeling of eigenspaces of the canonical momenta tTx 
and TTy by irreducible projective representations of T. Substituting P — —j-tt' ■ R one can write 

r(R)=expP = ^ 



n< 

n=0 



and 



The last commutator equals 



where 



OC _j 

[expP,7re] = ^-[P",^5], ^^x,y. (8) 



n 

n=0 



ri-1 



Cf = [P, J and = TT^ . 



These formulae yield that = [P, ttj] = is a sufhcient condition for [T(R),H] = 0. For example, it can be solved 
in the case of the homogeneous magnetic field, i.e. when tt is a linear function of x and y, and this solution was found 
by Brown and Zak and then generalized to any linear vector potential by the author (see also ) . However, the case 
of the inhomogeneous magnetic field yields A being a square (or higher order) function of coordinates x, y and the 
condition Cf = does not lead to any non-trivial solutions for P. Therefore, the formula (^ has to be substituted 
to the condition (^) and the solution has to be found in this more general case. So, one obtains 

[-P p- = E ^ E ( : ) cL.p' = f E >l) -p p • (10) 

n=0 ■ fe=0 ^ ^ \n=l ■ / 

Therefore, the sum of commutators C| has to be equal to zero, what, is satisfied by the special solution — 0. After 
substitution tt^ — —ihd^ — {e/c)A'^ the operator P = —j-Xn'^ — j-YiTy, where X = rixa, Y = UyU are coordinates of 
R, decomposes into two summands: 

ig 

P^Pi+P2; Pi^^{Xdx + Ydy), P2^—{XA'^ + YA'y). (11) 
So, in a general case, we have 



Cl = [P, TT,.] = - [{Xd, + Ydy)A^ - d^iXA', + YA'y)] . (12) 



c 



If one assumes that dxA^ = dyAy — 0, this can be simplified to 

Cf = -[YdyAx^dx{XA'^ + YA'y)] 
c " 



Cf = -^[XdxAy-dy{XA'^ + YA'y)]. (13) 



Example Let us consider the case A^ = ay and Ay — f3x + Hix /2. We are looking for a solution in a similar 
form, i.e. determined by an associated vector potential A' being a square function of x and y, so we assume that 

The equation ( p^ ) gives 



A'^ = a'y + a"y'^ and A'y = p'x + (i"x'^ 



Cf = -Y{a -(3' - 2(5" x) ; 
c 

Cf = -^X{f3 + Hix-a' -2a"y). 
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Since 7^ in a general case, we have to calculate C| (note that only Pi is relevant): 

Q = [P,C^] = 2P"-XY; 
C« - [P, Cf ] = ^X{-XHi + 2Ya") . 

Therefore Ci = for / > 3 and the commutator [expP, tt^] = if + C|/2 = for any x^y,X,Y. For ^ = a; we 
obtain 

/?" = and 13' ^ a. 

It is easy to notice that a solution for £, = y and Hi ^ Q can be obtained only when a" = and X = Q. It means 
that (magnetic) translations r([0, K]) commute with the Hamiltonian. It is not surprising since the magnetic field 
H = {[3 — a) + Hix is not periodic in the x-th direction. 

This example shows that the linear term in A^^ (the constant term in H) always appear in A'^ in the same way 
as in the case of the homogeneous magnetic field. Therefore, for the vector potential given by Eq. (||) we choose the 
associated potential in the following form 

A'^ = py + f{x,y), 

A'y = ax + g{x,y), (14) 
where functions / and g will be determined from the commutation conditions. Substituting it to Eq. (|l^) we obtain 

Cf = -^ {YHiCos{ky)+d.AXfix,y) + Ygix,y)]} ; 

Cf = ^{XHiCOs{kx)~dy[Xfix,y) + Yg{x,y)]} . (15) 



The next commutators will determined with the use of the operator Pi , see Eq. ( |11[ ) , and derivatives of each function 
will be consider separately. The first part of Cf , up to the constant factor —Hie/kc, gives the following series 

kYcos{ky), (kY)^ sin{ky) , -{kYf cosiky) , -(fcr)^ sin(fcy) , {kYfcos{ky),... 

Decomposing it to two series with cos{ky) and sm{ky), respectively, and substituting it to the infinite sum in Eq. (^^ 
one obtains (taking into account that k — 2t: /a and Y — Uya, so kY = 2nj,7r) 

, (-i)"(/fcr)2"+i ^, / ^ (-i)"(A:y)2»\ 

= cos(fcy) sin(fcr) + sin{ky){l - cos{kY)) ^ . (16) 

The same result is obtained for a part of Cf containing cos(fca;). Therefore, operators T(R) commute with the 
Hamiltonian for trivial functions 

f{x,y) = 9{x,y) = 0. 

It is not so surprising if we recall that the factor system depends on a magnetic flux through a lattice cell [l5[| and 
the periodic part of H gives no input to it. Such solution leads to interesting difference between H and H': the first 
is periodic, whereas the second is uniform since H' = —Hq. This difference is more evident if we take into account 
the vector potentials A and A': the original one satisfy an inhomogeneous wave equation 

4:71 

V^A = J = kHi[sin{ky), - sin[A;a;)] , 

c 

whereas V^A' — 0. 

The obtained solution means that a spatially periodic magnetic field does not require any changes in the group 
of magnetic translations (projective representations of the translation group). 
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3.4 Other solutions 



It is interesting to check what happen in the case when A'^ depends on x and ?; in a similar way as A^. To begin 
with we choose 

f{y) = -y sin(fcy) , g{x) = — sm{kx) . (17) 
Such a choice leads to the associated magnetic field H' equal to 

H' = {a~ P)- Hi [cos{kx) + cos{ky)] = -H , 

as we have obtained in the case of the homogeneous magnetic field (c/. fl^). The series of commutators Cf containing 
g{x) and its derivatives is calculated as follows (up to the constant factor HiYe/kXc) 

kXcos{kx), (kXysinikx), ~{kXfcos{kx) , -(kX)^ sm{kx) , {kXfcos{kx), ... 

Substituting it to the infinite sum in Eq. ( |l0| ) one obtains, as in Eq. (|l6|), 

cos(fcy) sin(fcy) + sin(fcj/)(l — cos(fcF)) — . 

However, this solution has been obtained after division by X, so it is valid only for X ^ 0. Analogous considerations 
for f{y) lead to the same result, but now Y ^ 0. Therefore, the potential obtained for functions ( pT| ) determines 
operators commuting with the Hamiltonian only ioi X,Y ^ or X = Y = Q. Exclusion of the axis X = is caused 
by the lack of one power of X. It can be revealed if the magnetic translations will be defined locally [ p^ , i.e. a form 
of functions / and g will depend on X and Y . The simplest solution is to put 

fY{y)=Yf{y) and gx{x)^Xg{x). 

Another way is realized by switching / and g, so 

H H 

f{x) = — sin(fca:) , g{y) = — ^ sin(fcy) . (18) 

In this case, however, H' = —Hq (without the periodic term) and, moreover the vector potential A' is no longer 
written in the radiation gauge. 

3.5 Movement of a Bloch electron 

The approximated considerations of an electron in spatially inhomogeneous field presented in many textbooks, e.g. 
|]l7| , show that for small values of \7H a charged particle gains an additional velocity perpendicular to H and VH. 
Since VH is perpendicular to lines on which H is constant then it is most likely that the orbit center will move along 
these lines. Results of simple numerical simulations (only the Lorentz force for a periodic magnetic field has been 
taken into accout ) seem to confirm this picture. However, a system is very unstable and actual behavior strongly 
depends on starting values and accuracy. When a time-step has been relatively large {uiAt ~ 0.01), the center moved 
towards the maximum (for large H) or minimum (for small H) of the function H(x, y). This movement was almost 
negligible for small values of At (cjAt ~ 0.0001). A critical value of H depends on a particle velocity and corresponds 
to the orbit diameter 2p close to the field period, cf. It could be expected that in the limit case A< At the 
orbit center would be 'stabilized' on an 'equipotential' line H —const. The other conclusions are almost obvious: 
(i) in a strong magnetic field (a very small orbit diameter) the movement of the orbit center is much slower since a 
particle moves in almost constant magnetic field; (ii) in the case of small magnetic field (a large orbit diameter) a 
particle moves in regions of different magnetic field magnitudes and the orbit center trajectory is not a smooth line 
(see Fig. 0). 

4 Final remarks 

It has been shown that a spatially periodic magnetic field leads to the same magnetic translation operators as in the 
case of a uniform field and a periodic potential V{y). Since the results obtained do not depend on an actual value of 
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Figure 1: A trajectory of the orbit center for a large orbit radius p ^ a 



k then they are vahd for any periodic function of x and y written as its Fourier transform. When a periodic potential 
V{r) is also present than our considerations are applicable to the case of commensurate periods of H and V only. 
The results of numerical simulations suggest that the orbit center moves along an 'equipotential' line H = const. This 
movement is more stable in the case of relatively large fields, i.e. for the orbit radius p smaller than the field period 
a. The magnetic translation operators cannot be interpreted as the position of the orbit center |]l5[ |l^, since the 
later is not a constant of motion. They rather correspond to the center of the orbit center trajectory. It should 
be underline that this work does not present a theory of the observed effects but provides a mathematical tool, 
which can be useful in such investiagtions. 
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